An Improved Model for Polyether Production from 1,3‐Propanediol

A dynamic mathematical model is developed for production of Cerenol polyether from 1,3‐propanediol in a batch reactor system. The model accounts for polycondensation reactions and side reactions in the liquid phase and for mass transfer of volatile species to the vapor. Parameters are estimated using measured liquid‐phase concentrations of monomer, oligomers, water, and end groups as well as the mass and composition of condensate collected from the overhead condenser system. The proposed model uses novel probability factors to keep the model equations relatively simple while accounting for the complex influence of superacid catalyst on reaction rates. The model is a significant advance over previous Cerenol models because it better accounts for mass‐transfer rates and for the dynamic behavior of the condenser. In addition, the proposed model accounts for the inhibitory influence of water on polycondensation kinetics due to hydration of hydroxyl ends. The model equations and parameter estimates provide a substantial improvement in fit to the data, especially for long reaction times and high catalyst levels, resulting in a 97% reduction in the value of the weighted least squared objective function compared to equations and parameters from a previous model.     

Stochastic method to control Mycobacterium tuberculosis epidemic

We study the origin of TB (tuberculosis) epidemic and complex distributions of various populations of TB infection within the stochastic framework. The stochastic nature of this disease infection could be linked to the stochastic behaviour at genome level which is exhibited in SNP (single nucleotide polymorphism) distributions of experimentally identified hotspot driver genes. Our results show the emergence of random clusters, and well-defined discrete domains of the respective species populations in the model driven by demographic stochasticity and intrinsic complex species interaction. The multifractal analysis of the time series of the species populations indicate that TB epidemic could be mainly caused by contact communication and is directional. We propose that any TB epidemic may have high chance of approximately periodic recurrence and can be controlled by optimizing some of the parameters involved in the system modelling.     

Numerical approach for the evolution of spin-boson systems and its application to the Buck–Sukumar model

We study the evolution properties of spin-boson systems by a systematic numerical iteration approach, which performs well in the whole coupling regime. This approach evaluates a set of coefficients in the formal expression of the time-dependent Schr?dinger equation by expanding the initial state in Fock space. This set of coefficients is unique for the spin-boson Hamiltonian studied, allowing one to calculate the time evolution from different initial states. To complement our numerical calculations, we apply the method to the Buck–Sukumar model. We find that when the ground-state energy of the model is unbounded and no ground state exists in a certain parameter space, the time evolution of the physical quantities is naturally unstable.     

Single top partner production in the tZ channel at eγ collision in the littlest Higgs model with T-parity

Introducing the top partner is a common way to cancel the largest quadratically divergent contribution to the Higgs mass induced by the top quark. In this work, we study single top partner production in the tZ channel at eγ collision in the littlest Higgs model with T-parity(LHT). Since it is well known that polarized beams can enhance the cross section, we analyze the signal via polarized electron beams,and photon beams. we have selected two decay modes for comparison, based on the leptonic or the hadronic decays of the W and Z from the top partner. We then construct a detailed detector simulation, and choose a set of cuts to enhance signal significance. For mode A(B), the capacity for exclusion in this process at s~(1/2)=3TeV is comparable to the current experimental limits with L=1000(500) fb~(-1). If the integrated luminosity can be increased to 3000 fb~(-1), the top partner mass+mTcan be excluded up to 1350(1440) GeV at 2σ level. We also considered the initial state radiation effect, and find that this effect reduces the excluding ability of the eγ collision on the the top partner mass by approximately 10 GeV. Moreover, the ability to exclude the LHT parameter space at eγ collision complements the existing research.     

Stochastic delay foraging arena predator–prey system with Markov switching

Abstract

In this paper, we introduce white noise, telegraph noise and time delay to the two-dimensional foraging arena population system describing the prey and predator abundance. The aim is to find out how the interactions between white noise, telegraph noise and time delay affect the dynamics of the population system. Firstly, the existence of a global positive solution is verified. Then the long-time properties including the stochastically ultimate boundedness, extinction and some other asymptotic pathwise estimation of this population system are studied. Finally, the main results are illustrated by two examples.     

具有非线性传染率和时滞的随机SIQR计算机病毒模型的渐近行为

研究了一类具有非线性发生率和时滞的随机SIQR计算机病毒模型.首先证明了该系统具有唯一的全局正解,然后通过构造适当的Lyapunov函数并利用伊藤公式,分析了该模型的解在无病平衡点附近及地方病平衡点附近的渐近行为,最后通过数值模拟对随机系统解的渐近行为做了进一步的分析并给出了结论.     

Readjustment of the Bohr stopping force from energies of keV/n to a few tens of MeV/n ions in elemental targets

An analytical approach for the electronic stopping force for non-relativistic energies that has no adjustable parameters has been developed. The approach combines the Bohr model for the close collisions and the Firsov model for the distant collisions. In order to combine the two models, a probabilistic model was introduced. We have applied our model to ¹⁶O in ¹²C, ¹⁶O in ²⁷Al, ⁸⁴ Kr in ²⁷Al, ⁵Li in ¹²C, ¹²C in ¹²C, and ¹³²Xe in ¹²C systems and compared with SRIM/MSTAR software, the original Bohr model, the Firsov model and available experimental data. We have found that the calculated electronic stopping force values are in agreement with the general qualitative behaviour of the electronic stopping force as a function of particle velocity reported in the literature. The proposed analytical formula is expected to be valid for other projectile-target combinations but more experimental data are needed to verify this assumption.     

基于激光雷达的垂直能见度反演算法及其误差评估

针对大气垂直方向上消光系数分布不均匀难以用传统方法直接测量垂直能见度的问题,提出了一种基于激光雷达探测垂直能见度的计算方法。根据大气辐射传输基本原理,借助于辐射传输方程,推导出了垂直能见度的计算公式;然后利用激光雷达原理方程和Klett算法反演出大气垂直方向上的消光系数分布,基于此提出了垂直能见度的迭代算法。最后,利用灰色模型GM(1,1)和批统计算法,对激光雷达反演得到的后向散射系数进行了评估,给出了误差置信区间为(0.760±0.339)×10^-4(srad·km)^-1。结果表明,该方法是一种特别有效的计算垂直能见度的方法,符合探测的基本需求,且误差小精度高。     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

带有潜伏期的狂犬病模型及控制策略研究

研究建立了狗与人两个群体的SEIR狂犬病模型,把狂犬病潜伏期阶段的传染性加入到模型中,利用遗传算法对1996年到2015年的数据进行拟合,估出了中国狂犬病的基本再生数R0≈1.66,证明了狂犬病无病平衡点的全局稳定性以及模型的一致持续性.根据数值模拟结果,发现处于潜伏期的犬类对人类狂犬病造成威胁.因此控制狂犬病的一个关键措施就是要做好犬的管理,减少与犬只的接触,这对抑制狂犬病的爆发有着重要的作用,也为我国狂犬病的控制提供了可靠的理论依据.